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■ Telescope. Since the NB816 imaging is so deep we can detect 10 well-defined NB816- 

Ph 1 dropout galaxies that we identify to lie at z — 5.8 - 6.5. We discuss their observational 
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^ ' 1. Introduction 



In order to search for high-redshift galaxies around z ~ 6 in the Subaru Deep 
Field, we have investigated NB816-dropout galaxies where NB816 is the narrowband 



Surveys for high-redshift galaxies provide observational bases to explore the formation 
and evolution of galaxies. A number of such surveys have been made in this decade. Indeed, the 
optical broad-band color selection technique has been utilized to find large samples of high- z 
galaxies at z ~ 3 - 7; i.e., Lyman break galaxies (LBGs) (e.g., Lanzetta et al. 1996; Madau et 
al. 1996; Steidel et al. 1999; Ouchi et al. 2004; Giavalisco et al. 2004; Stanway et al. 2004; 
Bouwens et al. 2004). Deep optical surveys with a narrowband filter have also been conducted 
to search for strong Lyman a emitters (LAEs) at high redshift. Since the early success of 
narrowband filter surveys (e.g., Hu et al. 1996; Cowie & Hu 1998), a large number of high- z 
LAEs have been found (e.g., Hu et al. 2002, 2004; Kodaira et al. 2003; Malhotra & Rhoads 
2001; Taniguchi et al. 2005; see for reviews, Taniguchi et al. 2003; Spinrad 2004). Although 
such deep narrowband imaging data are basically used to find strong emission-line objects, they 
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can also be used as one of the photometric data points. For example, Cuby et al. (2003) made 
their optical imaging survey for LAEs at z ~ 6.6 by using both broad-band data (B, V, R, 
and /) and narrow-band data taken with their NB920 filter. They found a LAE at z = 6.17 
from a sample of NB920-excess objects. In this case, the NB920 data were used to sample the 
continuum above the Lyman break. This successful result suggested that any narrowband data 
can be used to find not only strong LAEs but also so-called dropout high- ,2 galaxies. 

Recently, Shioya et al. (2005) proposed a new method, intermediate-band dropout 
method, to search for high- 2 galaxies (see also Fujita 2003). In this method, an intermediate- 
band filter at A > 7000 A is used together with typical i' and z' filters. We focus on dropouts 
using the intermediate-band filter, while Cuby et. al. (2003) selected a z=6.17 Lyman break 
galaxy using a narrowband filter. The merit of the new method is to distinguish either very 
late-type stars such as L and T dwarfs or dusty galaxies at intermediate redshift from real high- z 
LBGs. The reason for this is that such interlopers do not show strong narrow-band depression 
although they have very red colors that are indicative of LBGs. In the survey method, we may 
use a narrowband filter at A > 7000 A (Hu et al. 2004). We also note that Y.Kakazu (private 
communication) proposed the NB816- depressed method to separate high-reshift quasars and 
galaxies from L and T dwarfs based on the same principle. In order to find high- z galaxies 
around z ~ 6, we use the data archive of the Subaru Deep Field (Kashikawa et al. 2004). Since 
in the SDF, deep imaging data with the narrowband filter NB816 are available, they are very 
useful to demonstrate that the narrowband dropout technique is powerful to find such high- z 
galaxies. 

2. The Narrow-Band Dropout Method 

Our new method is basically to find narrow-band-dropout objects from a deep imaging 
survey data set. In order to describe our method practically, we consider the case of NB816 
dropouts where NB816 is a narrow-band filter centered at A c ~ 815 nm with a FWHM = 12 
nm (e.g., Taniguchi et al. 2003b and references therein; Hu et al. 2004). 

In this case, we need deep imaging data taken with filters of i', z', and NB816. In Fig. 
1, we show the transmission curve of the three filters taken from the data for Suprime-Cam 
filters 1 (Miyazaki et al. 2002). 

We also show spectral energy distributions (SEDs) of a LBG (starburst) at z ~ 5.9, an 
L dwarf star in our Galaxy, and a dusty galaxy at z ~ 1.1. The SED of a LBG is produced 
by GALAXEV (Bruzual & Chariot 2003). The SEDs of local galaxies are well reproduces 
by models whose star-formation rate declines exponentially (the r model): i.e., SFR(t) oc 
exp(— t/r), where t is the age of galaxy and r is the time scale of star formation. We set r = 1 
Gyr and t — 1 Gyr to derive the SED of a starburst (SB) galaxy. To make an SED of an SB 

1 http://www.subarutclescope.org/Observing/Instrumcnts/SCam/scnsitivity.htnil 
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with emission lines, we calculate emission-line fluxes using the formula of L(H/?)(ergs s _1 ) = 
4.76 x 10 _13 A r L yc , where N^ yc is the ionizing photon production rate (Leitherer & Heckman 
1995), and the table of relative luminosity to H/3 luminosity in PEGASE (Fioc & Rocca- 
Volmerange 1997). The SED of an LBG is the SB model with Ay = and that of a dusty 
galaxy is the SB model with Ay = 10 mag adopting the reddening curve of Calzetti et al. 
(2000). It is shown that the LBG at z = 6 is selected as a NB816 dropout. However, the L 
dwarf is not selected as a NB816 dropout even though it shows a very red %' — z' color. In this 
way, we can select LBGs at z > 5.8 as NB816 dropouts. 

In Fig. 2, we show the diagram between i 1 — z' and NB816 — z' for galaxies found in 
the Subaru Deep Field (Kashikawa et al. 2004). The selection procedure is given in the next 
section. Galactic stars (Gunn & Stryker 1983) and L and T dwarfs 2 (Kirkpatrick 2003) are 
also shown. An evolutionary path of a LBG is shown by the blue curve; its SED is generated 
by using the GALAXEV (Bruzual & Chariot 2003). 

When we select LBG candidates at z ~ 6, we adopt a color criterion of i' — z' > 1.3 - 1.5 
(e.g., Nagao et al. 2004; Stanway et al. 2004; Giavalisco et al. 2004). However, this criterion 
is not enough to select such LBGs unambiguously because L and T dwarfs also have such very 
red colors. In addition, some dusty galaxies at intermediate redshift may also have such red 
colors (e.g., a part of extremely red objects; see Hu & Ridgway 1994). 

Here we would like to stress that LBGs at z > 5.8 are efficiently selected if we also use 
another criterion of iV.B816 — z > 1.5-2; i.e., NB816-depressed objects. The most important 
merit of this method is that LBG candidates are selected more unambiguously because this 
new criterion can separate LBGs from interlopers (i.e., L, T dwarfs and dusty galaxies at lower 
redshift). Hu et al. (2004) noted that LAEs at z ~ 5.7 are well separated from L and T dwarfs 
if a proper criterion on the NB816 — z' color is adopted. 

3. Application to the Subaru Deep Field Data 

In order to demonstrate how the new NB dropout method works well, we apply it to the 
Subaru Deep Field data (Kashikawa et al. 2004). Photometric catalogs (Version 1) have been 
made public 3 . In the following analysis, we use a 2"0 aperture magnitude and the error in 2"cf) 
aperture magnitudes from the catalog. We also calculate errors of colors using the errors in the 
2"(f) aperture magnitude. 

As shown in Fig. 2, being independent of the strength of the Lyman a emission, galaxies 
at z > 5.73 can be identified as NB816-dropout objects. Assuming that there are no elliptical- 
like red galaxies at z ~ 2, we estimate a reddest NB816 — z' color of ~ 1.3, corresponding to 
the color of elliptical galaxies at z ~ 1. Therefore, we adopt the following selection criterion for 

2 http :/ /spider . ipac . caltech . edu/st aff / davy / ARCHIVE / 

3 http://soaps.naoj.org/sdf/ 
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NB816-dropout galaxies, 

NB816-z' >1.5 (1) 

and 

zf< 26.07 (5a). (2) 

In order to reduce any contaminations from foreground objects that are free from ab- 
sorption by the intergalactic neutral hydrogen, we also adopt following criteria 

B> 28.45 (3a), (3) 
V > 27.74 (3a), (4) 
and i?> 27.80(3(7). (5) 

By using the above three criteria, we first obtain a sample of 42 NB816-dropout objects. From 
42 NB816-dropout objects, we select the objects with iVi?816 — z' > 1.5 above the 3a error as 
the more reliable NB-dropout objects. The number of more reliable NB816-dropout objects is 
15. 

Next, we must separate high- z galaxies from Galactic low-temperature stars (L and T 
dwarfs). It is known that most L and T dwarfs satisfy the following relation: 

(i 1 - z') - (NB816 - z') ~ 0.8, (6) 

as shown in Fig. 2. Taking account of the scatter of the colors of L and T dwarfs, we can 
separate LBG candidates from them by using the following criterion, 

(if -z')- (NB816-Z 1 ) < 0.5. (7) 

In order to securely select a reliable sample of NB816 dropout galaxies, we pick up objects that 
satisfy the condition (7) above the 3a error conditions from 15 more reliable NB816-dropout 
objects. Then we finally obtain a sample of 10 NB816-dropout galaxies. Their photometric 
properties are summarized in Table 1. In this table, we present 3a error of colors. 

We note that galaxies with faint continuum and strong-emission lines in z', e.g., strong 
[OIII] emission-line objects at z ~ 0.8, may satisfy the above selection criteria. Actually, Y. 
Kakazu has found strong [OIII] emission-line objects in her NB816-depressed sample (private 
communication). Since there is an NB921 observation for SDF, we can separate such objects 
from genuine NB816 dropouts using z' — NB921 colors. If strong-emission lines lie in both z' 
and NB921, those galaxies would be observed as iV.B921-excess objects. If strong-emission 
lines do not lie in NB921 but in z', iV\B921-band of those galaxies must be very faint. All 
our NB81 6-dropout sample are detected (> 3a) in NB921-band, and their zf — NB921 colors 
are nearly flat, which range from —0.39 to 0.02. Therefore, they are considered to be galaxies 
at redshift between 5.73 and 6.50. One object, ID 35247, is considered to be a LAE since 
its z! - NB921 color is 1.49. In fact, this object is identified as SDF J132408.3+271543 in 
Taniguchi et al. (2005) and confirmed as a LAE at z — 6.554. 



Now, we consider the UV luminosity function of star-forming galaxies at z ~ 6. In 
addition to the UV luminosity functions of star- forming galaxies at z ~ 3 to ~ 5 (e.g., Steidel et 
al. 1999; Ouchi et al. 2004), those of star-forming galaxies at z ~ 6 were also recently obtained 
by red %' — z' color (i'-dropout galaxies, e.g., Bouwens et al. 2004) or by using the %' — Zr versus 
zb — zr diagram (Shimasaku et al. 2004). Hu et al. (2004) also obtained the UV luminosity 
function of Lyman a emitters (LAEs) at z ~ 5.7. 

Since the observed ^'-magnitude of a galaxy depends on not only the UV luminosity 
but also the redshift of the galaxy and the strength of the Lya emission line, we compare our 
results with the prediction evaluated under some assumptions. If the UV luminosity function 
of galaxies at z ~ 6 is the same as that at z ~ 4.0 or 4.7 (Ouchi et al. 2004), we can evaluate 
the number of star-forming galaxies at z ~ 5.8 to 6.6. We show the expected galaxy number 
evaluated using the SB model SED without emission lines in Fig. 3. We also show the expected 
number of galaxies for the UV luminosity function of LBGs at z ~ 6 (Bouwens et al. 2004) and 
that of LAEs at z ~ 5.7 (Hu et al. 2004). Since the strong Lya emission line modify the z' 
magnitude for galaxies at redshift higher than « 5.75, we use our SB model with emission lines 
for the UV luminosity function of LAEs. The rest-frame equivalent width of Lya emission is 
about 70 A in our SB model with emission lines. Taking account of the completeness for the 
SDF sample as a function of z' magnitude (Kashikawa et al. 2004), we expect that the number 
of galaxies brighter than z' = 26.0 ranges from 42 to 200 (Fig. 3). 

The selection criteria of our final sample are so severe that the number of our final sample 
(10) may give a lower limit. On the other hand, the number of our first sample (42) may give 
a upper limit, since it may include contaminations. Fig. 3 suggests that the number of LBGs 
at z ~ 6 is smaller than that at z ~ 4 - 5. This trend is consistent with that derived from the 
observation of the Hubble Ultra Deep Field (Bunker et al. 2004) and the Hubble Ultra Deep 
parallel fields (UDF PFs: Bouwens et al. 2004). Although the reason why the number of our 
sample is smaller than the prediction using the UV luminosity function at z ~ 6 is still not 
clear, we note that the number density of bright LBGs at z ~ 6 in the SDF is 60 % of the value 
obtained in the UDF PFs by Bouwens et al. (2004) (Shimasaku et al. 2005). 

4. REMARKS 

In Section 3, we have demonstrated that the NB816 dropout method is capable of 
selecting reliable candidates of LBGs at z > 5.8. Although we have adopted the NB816 dropout 
method there, our method is also useful when we use other narrow-band filters. 

For example, NB921 is useful to find galaxies at z > 6.7 while NB711 is useful for galaxies 
at z > 4.9. It is also noted that a rest-frame Ha filter (NB656) and a rest-frame [O ill] filter can 
be used to search for galaxies at z > 4.4 and at z > 3.1, respectively. Since there are a number 
of deep survey data in which some of narrow-band filters are used, it is recommended that the 
narrow-band dropout method is applied to them, providing new samples of LBGs/LAEs at 
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high redshift. 
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Fig. 1. The upper panel shows the transmission curves of the filters, i', z', and NB816. The lower panel 
shows SED of Lyman break galaxy (LBG) at z ~ 5.9 (blue), dusty starburst at z ~ 1.1 (red), and Galactic 
L dwarf stars, 2MASS J09510549+3558021 (green). The SEDs of LBG and dusty starburst are produced 
by GALAXEV (Bruzual, Chariot 2003, see text). Although, these objects shows the same i' — z' colors 
(~ 2). The NB816 — z' colors of LBG is redder than those of a L dwarf star and a dusty starburst galaxy 
(see Fig. 2). 
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Fig. 2. i' — z' vs. ^58168 — z' diagram displaying the colors of model galaxies and stars. Solid lines arc 
loci of model galaxies produced by GALAXEV (Bruzual, Chariot 2003): red line - r = 1 Gyr model with 
age of 8 Gyr, green line - r = 1 Gyr model with age of 1 Gyr (SB model), and blue line - SB model with 
emission lines (see text). Orange marks show stars. Asterisks are Gunn & Stryker (1983)'s star, crosses 
are L dwarfs and pluses are T dwarfs. Open squares with error bars show our final sample. 
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Fig. 3. Number of NB816-dropout galaxies as a function of z' magnitude. Triangles show the number of 
our final sample and inverse triangles show the number of first sample. Lines show the expected number of 
galaxies with z > 5.73 for different assumptions. Thin dotted (dot-dashed) line shows the case for the SED 
without Lya emission together with the UV luminosity function at z ^4.0 (4.7) (Ouchi ct al. 2004). Thick 
solid line shows the case for the SED without Lya emission together with the UV luminosity function 
at z ~ 6.0 (Bouwens et al. 2004) and thick dashed line shows the case for the SED with Lya emission 
together with the UV luminosity function at z ~ 5.7 (Hu et al. 2004). 
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Table 1. Photometric properties of NB816-dropout galaxies 



ID 


Name 


Optical AB magnitude 


Optical colors 






i' z' NB816 NB921 


i'-z' NB816-z' 



18309 


SDF 


J132436 


.6- 


-271246 


27, 


93 


25 


.89 


29 


.67 


26 


.01 


2 


nc+0.47 
UO -0.33 


3 


■78±?°75 


35247 


SDF 


J132408 


.3- 


-271543 


27. 


.62 


25 


.98 


28 


.37 


24 


.49 


1 


■ o4 -0.27 


2 


■39ta78 


39694 


SDF 


J132418. 


,4- 


-271633 


27 


33 


25 


.77 


27 


,99 


25 


.75 


1 


C13+0.28 
' OO -0.22 


2 


99 + 1.47 
' zz -0.60 


46427 


SDF 


J132444. 


,8+271749 


27 


,82 


25 


.95 


28 


.20 


25 


.96 


1 


■ 8 '-0.31 


2 


z ^-0.69 


66699 


SDF 


J132406. 


,0+272222 


27 


,86 


25 


.90 


28 


.79 


26 


.00 


1 


qfi +0.44 
aD -0.31 


2 


.90t 1 To2 


77541 


SDF 


J132344. 


.7- 


-272450 


27 


.52 


26 


.05 


28 


,84 


26 


.06 


1 


zlo+0.35 
■ 48 -0.26 


2 




94464 


SDF 


J132452. 


.1- 


-272842 


27 


.34 


25 


.34 


28 


,27 


25 


.23 


2 


no+0.26 
uz -0.21 


2 


no+3.29 
' JO -0.73 


106928 


SDF 


J132339 


.7- 


-273138 


27 


.44 


25 


.49 


28 


.12 


25 


.66 


1 


yD -0.23 


2 


fi o + 1.94 
°°-0.66 


117034 


SDF 


J132436 


,0- 


-273356 


26 


,50 


25 


.70 


28 


,40 


25 


.93 





on+o.i8 


2 


70+°° 
■ ' u -0.79 


176181 


SDF 


J132507 


.7- 


-274509 


28, 


.16 


25 


.28 


28, 


,54 


25 


.49 


2 


on+0.56 
oy -0.37 


3 


97+00 
■ z '-0.87 
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